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> ■ ABSTRACT 

O ' 

#y^ ' PG 0014+067 is one of the most promising pulsating subdwarf B stars for seismic analysis, 

fT^ . as it has a rich pulsation spectrum. The richness of its pulsations, however, poses a fundamental 
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challenge to understanding the pulsations of these stars, as the mode density is too complex to be 
explained only with radial and nonradial low degree {(, < 3) p-modes without rotational splittings. 
One proposed solution, suggested by Brassard et al. (2001) for the case of PG 0014+067 in 
particular, assigns some modes with high degree (^ = 3). On the other hand, theoretical models 
of sdB stars suggest that they may retain rapidly rotating cores (Kawaler & Hostler 2005), 
and so the high mode density may result from the presence of a few rotationally-split triplet 
(£ = 1), quintuplet {t = 2) modes, along with radial {£ = 0) p-modes. To examine alternative 
theoretical models for these stars, we need better frequency resolution and denser longitude 
coverage. Therefore, we observed this star with the Whole Earth Telescope for two weeks in 
October 2004. In this paper we report the results of Whole Earth Telescope observations of the 
pulsating subdwarf B star PG 0014+067. We find that the frequencies seen in PG 0014+067 do 
not appear to fit any theoretical model currently available; however, we find a simple empirical 
relation that is able to match all of the well-determined frequencies in this star. 

Subject headings: stars: subdwarfs — stars: oscillations — stars: individual: PG 0014+067 



1. Introduction 

Even though stellar evolution is one of the most 
mature areas in theoretical astrophysics, there are 
still phases in the life of an ordinary star (like our 
Sun) that are not well understood. We have a 
nearly continuous evolutionary path of a low mas- 
sive star all the way from its birthplace, a proto- 
stellar cloud, to its graveyard, with few gaps in 
our understanding. However, at least one major 
gap remains: the core helium flash that initiates 
helium core burning. Unlike most phases of stel- 
lar evolution, the flash itself occurs in a matter 
of minutes rather than millions of years. It is es- 
sentially a massive runaway nuclear reaction near 
the center of the star that releases a huge amount 
of energy on a short time scale. Despite ignition 
on a dynamical time scale, apparently this event 
does not make the star explode, because we do 
see stars that are clearly survivors of this process. 
Post helium core flash objects are identified with 
horizontal branch stars in globular clusters, and 
their field counterparts. 

At the extreme blue end of the horizontal 
branch, subdwarf B (sdB) stars are survivors of 
the core helium fiash that have a small hydro- 
gen rich envelope surrounding the helium-burning 
core. According to evolutionary calculations, for 
example those by Dorman ct al. (1993), sdB stars 
are core helium burning stars with a canonical 
mass of M w 0.5 Mq and a very thin, inert hydro- 
gen envelope with ^f^- ~ 0.0004 which places them 
on the hot end of the horizontal branch (HB), the 
so-called extreme horizontal branch (EHB). The 
fact that sdB stars have lost almost all of their 



hydrogen layer is what makes them so enigmatic 
from the stellar evolution point of view. To lose 
that much mass, they have to suffer a consider- 
able mass loss during the RGB and most probably 
during the helium core flash. 

Details of the origin of the sdB stars are largely 
unknown. There are currently many competing 
theories for their origin which involve several pos- 
sible channels. These include single star evolution 
with an extreme mass loss on the RGB (D'Cruz et 
al. 1996), or different binary interactions: common 
envelope binary evolution (Mengel et al. 1976), 
stable Roche lobe overflow (Han et al. (2002, 2003) 
and references therein) , or merging of two helium 
white dwarfs (Iben 1990). 

The most fundamental missing piece to our un- 
derstanding of the evolution of the sdB stars is the 
nature and physics behind their mass loss (Fusi- 
Pecci & Renzini 1976). To end up on the EHB 
they must lose nearly all of their hydrogen at al- 
most exactly the same phase as the helium core 
has attained the minimum mass required for the 
helium fiash to occur. By exposing the interior 
properties of these survivors we can learn how the 
stars manage to survive this event of unstable he- 
lium ignition and plug this gap in the story of stel- 
lar evolution. 

Luckily, some of these stars pulsate, letting us 
therefore use the tools and techniques of astero- 
seismology to probe their interiors. Discovered by 
Kilkenny et al. (1997), the pulsating sdB stars (of- 
ficially designated V361 Hya stars, but also called 
sdBV stars by analogy with the pulsating white 
dwarfs) are short period, low amplitude multi- 



mode pulsators. Typical periods are about 100- 
250s with an overall period range of 80-600s. The 
pulsation amplitudes are generally less than a few 
hundredths of a magnitude. The short periods, 
being of the order of, and shorter, than the ra- 
dial fundamental mode for these stars, suggest 
that the observed modes are low-order, low-degree 
p— modes (Charpinet et al. 2000). 

The rapidly pulsating sdB stars occupy a region 
in the logg — Tes plane with effective temperatures 
of Tcff = 29 000 - 36 000 K and surface gravities 
log g = 5.2 - 6.5. Many of the high gravity pulsat- 
ing sdB stars show complex pulsations that require 
a large number of modes. Radial modes alone are 
not enough to explain the complex light curves of 
these stars, therefore nonradial modes must also 
be present (Kawaler 1999; Kilkenny et al. 1999; 
Charpinet et al. 2000). 

Among the sdBV stars PG 0014+067 stands 
out as one of the richest pulsators. It was found 
to be an sdBV star in June 1998 (Brassard et al. 
2001). While close binaries hke KPD 1930-^2752 
(Poindexter et al. 2003) and lower-gravity sdB pul- 
sators like PG 1605+072 (Kilkenny et al. 1999) can 
show dozens of modes, PG 0014+067 is a higher 
gravity single sdB star that shows a large num- 
ber of nonradial modes, with at least sixteen re- 
ported in the discovery paper. Periods range from 
80 to 170 seconds. It is a "typical" sdBV star 
with an effective temperature of Tcff- — 33 310K 
and \ogg — 5.8 (Brassard et al. 2001). The rich- 
ness of its pulsation, however, makes it a star of 
vital interest to the study of sdBV stars and, by 
extension, the entire hot horizontal branch phe- 
nomenon. The rich mode density in PG 0014+067 
poses a fundamental challenge to understanding 
the pulsations of these stars, as there are not 
enough radial and nonradial modes with low de- 
gree £ {£ < 3) to account for all of them (without 
rotational splitting). Low degree nonradial modes 
have a limited frequency distribution which im- 
plies that we need to resort to higher values of 
£ to account for all of the observed periodicities. 
That presents further challenges, as geometrical 
cancellation increases rapidly with £ so that higher 
£ modes need to have a much larger intrinsic am- 
plitude than low £ modes to reach the same ob- 
servable amplitiude. 

One proposed solution, by Brassard et al. 
(2001), presents an envelope structural model that 



approximates their determined pulsation frequen- 
cies with p-modes of degrees 0, 1, 2, and 3. Their 
resulting seismic model provides an estimate of the 
mass of PG 0014+067 of 0.49 ± 0.02 M© and, in 
addition, matches the observed spectroscopic pa- 
rameters, and places strict limits on the thickness 
of the surface hydrogen layer. 

Another possible explanation for the mode den- 
sity in PG 0014+067 comes from theoretical cal- 
culations of the evolution of stars on the red giant 
branch, and the subsequent structure on the hor- 
izontal branch. These models show that stars can 
develop a rapidly rotating core on the red giant 
branch, and then preserve that core on the hor- 
izontal branch. If that is the case, then certain 
modes of pulsation in sdBV stars could show large 
rotational splitting - much larger than suggested 
by the (slow) surface rotation velocity. Models by 
Kawaler & Hostler (2005), in particular point out 
that rotational splittings of low order £ ~ I and £ 
= 2 modes can cause sufficient frequency peaks in 
the observed range provided that the stellar core 
rotates much faster than its surface. These models 
show rotational splitting that can vary by factors 
of two or more from one mode to the next. If that 
is the case in PG 0014+067, then it could suf- 
ficiently complicate the observed mode spectrum 
with modes oi £—0, 1, and 2 alone. 

The only published observations of PG 0014+067 
are single-site data from Brassard et al. (2001), 
who have observed this star using the 3.6 meter 
Ganada-France-Hawaii Telescope (GFHT) during 
five consecutive nights (1.9h average run length) 
in June 1998. These data show that the oscil- 
lation spectrum is complex, but being single-site 
data they suffer from severe one cycle per day 
aliases making the frequency identification rather 
difficult. 

To fully exploit this star we need to determine, 
unambiguously, the oscillation frequencies present. 
This can be achieved only by extended longitude 
coverage with a sufficient frequency resolution to 
be able to resolve the real pulsations from the 
inevitable one cycle per day aliases. With a se- 
cure set of oscillation frequencies, we can criti- 
cally examine the proposed theoretical models for 
these stars without concern for fitting artifacts 
caused by diurnal aliasing instead of true oscilla- 
tion modes. PG 0014+067 was selected as the pri- 
mary target for a Whole Earth Telescope (WET) 



campaign, held in October 2004. In this paper, 
we present the results of this multisitc campaign. 
In §2 we present the observations followed by a 
review of the data reduction procedure. In Sec- 
tion 3, we report on the frequency analysis of the 
WET data set on PG 0014+067, and we discuss 
the results in §4 along with the plans for future 
work. 

2. Observations 

PG 0014+067 was a primary target during the 
24th WET run, Xcov24, held in 2004. Being a 
faint star, with V ~ 16.5, and a low amplitude 
pulsator - most modes have amplitudes below 2 
mmag - we sought larger than average (>2 me- 
ter class) telescopes distributed around the Earth. 
Time was allocated to five such telescopes (> 2m) 
that were well-distributed in longitude. Obser- 
vations spanned two weeks centered around New 
Moon from October 7-21, 2004. With a right as- 
cension near Oh, and a near-equatorial declination 
of +7 degrees, PG 0014+067 was equally accessi- 
ble from both Northern and Southern Hemisphere 
sites during this time of the year. As a contin- 
gency, smaller telescopes (< 2 m) with CCD de- 
tectors were included to minimize the gaps in the 
data and to simplify the alias pattern of our spec- 
tral window. Together with 13 smaller telescopes, 
18 sites participated in this campaign. 

Unluckily, non-photometric conditions all over 
the globe during the main campaign resulted in 
less-than-ideal coverage. Table 1 lists all the obser- 
vations gathered during Xcov24 on PG 0014+067. 
As can be noticed from the log of the observa- 
tions given in Tabic 1, much of the data did not 
come from our bigger telescopes (> 2m), because 
of the bad weather conditions, but from moderate 
size ones (between Im and 2m). The 1.0 ni tele- 
scope at Lulin Observatory in Taiwan contributed 
data on six nights. An additional four nights of 
high-quality data came from the Bohyunsan Op- 
tical Astronomy Observatory (BOAO) 1.8m tele- 
scope in South Korea. Also, surprisingly good 
data were obtained by the SARA 0.9m telescope 
at KPNO, given the relatively small size of the 
telescope. Despite the terrible weather and some 
unexpected instrumental problems we had during 
Xcov24, the two week core campaign resulted in 
about 180 hours of data on the primary target. 



yielding a duty cycle of about 58 %. The WET run 
on PG 0014+067 spans two weeks which results in 
a formal frequency resolution of Aiy — 0.9/iIIz. 

The telescopes used for the primary target 
range in aperture from 0.9m to 2.56m. Data from 
NAOC (China), Moletai Obs. (Lithuania) and 
Pic Du Midi (France) were obtained using three- 
channel PMTs as described by Kleinman ct al. 
(1996). Exposure times were 5 seconds for 2m tele- 
scopes (China and France) and 10 seconds for the 
1.65m telescope in Lithuania. Channel 1 measured 
the program star, channel 2 measured a local com- 
parison star and the third channel simultaneously 
recorded sky background. All PMT photometric 
observations were taken in white light to maximize 
the photon count rate. 

All other data were obtained using CCDs (see 
column Instrument of Table 1). For the CCD mea- 
surements, exposures ranged from 5-15 seconds de- 
pending on the size of the telescope. As the read- 
out times differ from chip to chip and were quite 
long for some sites, observations were made in the 
"windowed" mode, i.e. only a portion of the chip 
was read out to minimize the cycle time, making 
sure there were always at least two comparison 
stars in the frame. Observations from McDon- 
ald observatory and LNA were made in the frame 
transfer mode. The actual CCD frames at every 
individual site were never bigger than 5x5 arc 
minutes centered around the target. The num- 
ber of comparison stars depended on the size of 
this window, having from two to six comparison 
stars in the frame. In this way, consecutive data 
points were obtained in 5 to 25 s intervals, depend- 
ing on the instrument. Data from McDonald and 
KPNO used BG 40 filters and data from NOT used 
a W#92 filter. Together with the science frames, 
each night of observation included calibration im- 
ages (bias, dark and flat field frames). 

Preliminary, "on the fly" data reduction was 
done during the campaign at a central headquar- 
ters (at Iowa State University). Since the WET 
collaboration does not operate with a standard 
CCD photometry extraction package, raw CCD 
frames were pre-reduced (overscan, bias, dark, flat 
filed corrected) by the observers using their own 
software immediately after the run, and the first- 
look photometry was sent to the headquarters 
(HQ) in the standard WET/XQED format. In- 
coming data were further reduced at HQ using the 



Table 1 
Journal of observations of PG 0014+067 in October 2004. 



Run name 


Telescope [m] 


UT Date 


Start [UT] 


Duration [hrs] 


Instrument 


mdr275* 


MDM 1.3 


05 


08:55:15 


2.10 


CCD 


teide-10072004* 


Tenerife 0.8 


08 


02:36:40 


0.55 


CCD 


sara213 


SARA 0.9 


08 


02:47:59 


8.74 


CCD 


PG0014_08oct04 


Loiano 1.5 


08 


19:24:50 


5.50 


CCD 


sara214 


SARA 0.9 


09 


01:58:59 


9.20 


CCD 


PG0014.09oct04* 


Loiano 1.5 


09 


19:45:09 


1.44 


CCD 


PG0014.10oct04 


Loiano 1.5 


10 


19:09:20 


6.35 


CCD 


lulin-0001 


Lulin 1.0 


11 


11:10:35 


8.35 


CCD 


041012A 


NAOC 2.16 


12 


12:27:00 


7.50 


PMT 


koroa-OOl 


BOAO 1.8 


12 


13:02:10 


5.45 


CCD 


korea-002 


BOAO 1.8 


13 


10:19:20 


5.32 


CCD 


041013A* 


NAOC 2.16 


13 


11:27:22 


8.50 


PMT 


lit-01 


Moletai 1.56 


13 


21:49:20 


4.65 


PMT 


francc-01* 


Pic Du Midi 2.0 


13 


20:50:00 


2.85 


PMT 


041014A 


NAOC 2.16 


14 


11:23:30 


8.49 


PMT 


koroa-003 


BOAO 1.8 


14 


10:55:16 


7.27 


CCD 


lit-05 


Moletai 1.56 


14 


19:00:10 


7.37 


PMT 


korea-004 


BOAO 1.8 


15 


10:19:00 


7.75 


CCD 


A0934 


McDonald 2.1 


15 


06:26:07 


2.65 


CCD 


NOT-001 


NOT 2.6 


15 


19:55:23 


4.32 


CCD 


NOT-002 


NOT 2.6 


16 


00:19:52 


4.87 


CCD 


lulin-0008 


Lulin 1.0 


16 


12:34:26 


6.45 


CCD 


NOT-003 


NOT 2.6 


16 


19:43:22 


8.37 


CCD 


rkbl* 


LNA 1.6 


17 


02:25:04 


3.15 


CCD 


A0939 


McDonald 2.1 


17 


05:19:58 


1.04 


CCD 


A0940* 


McDonald 2.1 


17 


07:32:18 


1.32 


CCD 


lulin-0010 


Lulin 1.0 


17 


11:36:36 


7.51 


CCD 


majal0182004 


KPNO 2.1 


18 


09:06:40 


1.78 


CCD 


lulin-0012 


Lulin 1.0 


18 


10:50:02 


3.10 


CCD 


041018B 


NAOC 2.16 


18 


13:47:00 


5.88 


PMT 


majal0192004 


KPNO 2.1 


19 


05:19:40 


3.40 


CCD 


lulin-0015 


Lulin 1.0 


19 


10:43:11 


8.33 


CCD 


majal0202004 


KPNO 2.1 


20 


02:11:20 


3.18 


CCD 


majal0202004b 


KPNO 2.1 


20 


09:39:40 


0.93 


CCD 


lulin-0016 


Lulin 1.0 


20 


11:35:48 


7.35 


CCD 


hungary-007* 


Konkoly 1.0 


20 


18:09:56 


1.56 


CCD 



*not used because of low signal-to-noise 



WET software XQED (Riddle 2003). These pre- 
liminary reductions were used during the course 
of the run to allocate telescopes between the pri- 
mary target and secondary targets, and for pre- 
liminary analysis of the pulsation spectrum of 
PG 0014+067. For this paper, all the data were 
carefully re-reduced by the lead author from the 
original CCD frames. 

For the PMT data, reduction followed standard 
WET procedures. After the proper dead time cor- 
rection, the sky background (channel three) was 
subtracted both from the target and from the com- 
parison star point by point. When necessary, the 
variable star data were then divided by the com- 
parison star to correct for transparency variations. 
To correct for any low- frequency variations caused 
by seeing and differential color extinction, the re- 
sulting light curve was then fitted with a low or- 
der polynomial. Bad points, coming from either 
electronic noise, guiding errors or clouds, were re- 
moved by hand. Finally, the light curve was mean 
subtracted leaving only the fractional variations 
from the mean {ma — A///) and a barycentric 
correction to the exposure mid-times was applied 
giving the times of the data points projected to 
the barycenter of the solar system. 

For our final reduction, the raw CCD data 
were reduced using the IRAF CCDPROC package 
and the reduction comprised correction for bias, 
dark counts and flat-field. Further photometric 
measurements (aperture photometry) on these re- 
duced frames were made using the PHOT pack- 
age. After carefully examining the extracted aper- 
ture photometry for a set of different apertures 
for each observation, the aperture that gave the 
best signal-to-noise ratio was used. The apertures 
that typically gave the optimal signal-to-noise ra- 
tio in the resulting light curve had a radius of 
«1.75 times the FWHM, with the surrounding an- 
nulus to determine the local sky level. In this way 
each star in the frame has its own corresponding 
sky channel. The aperture photometry measure- 
ments derived in this way were then imported into 
XQED. The resulting light curves, showing a frac- 
tional variation of the intensity of light as a func- 
tion of time are presented in Figure 1. 



3. The frequencies of PG 0014+067 

Combining all the light curves described in 
the previous section (Figure 1) we computed the 
Fourier transform (FT) to deduce the periodicities 
present in the data. To optimize the number of 
detected intrinsic frequencies present in the total 
light curve, gathered with the different aperture 
telescopes (from 0.9m-2.56m), we weighted the 
data using the weighting scheme favored by Han- 
dler (2003). Weights are determined for individual 
runs, with the weight for a run taken as inversely 
proportional to the mean point-to-point scatter 
within the run. However, to some degree, weight- 
ing inflates the amplitude of the daily aliases and 
widens the peaks. Consequently, weighting may 
affect the ability to resolve closely spaced peaks. 
Bearing this in mind weighted data were used to 
expose and identify periodicities, but the frequen- 
cies, amplitudes and phases of the individual peaks 
were determined using the 'unweighted' data. The 
weighted FT of the full data set on PG 0014+067 
gathered during Xcov24 is depicted in Figure 2 
showing the 'modulation amplitude', in units of 
mma= ma/1000, where nia=A///, of variations 
in the detected intensity / at a given frequency. 

To see the effect weighting has on the amplitude 
spectrum we have calculated the spectral window 
for the weighted and unweighted data. Figure 3 
shows the spectral window of the whole WET run 
on PG 0014+067 in both cases. The central peak 
in the plot corresponds to the input frequency of a 
single sinusoid; all other peaks are aliases caused 
by the gaps in the data set. 

3.1. Frequencies found in the Xcov24 data 

As can be clearly seen in Figure 2, the main 
action is concentrated in a very narrow frequency 
range. It appears as if all the power is distributed 
into two to three modes around the 7000/iHz with 
maximum intensity variations of ^ 3 mma. This is 
a typical signature of a low amplitude sdB pulsator 
(Kilkenny 2002). Figure 2 confirms that there are 
no frequencies with amplitudes above the mean 
noise level at higher frequencies, up to the Nyquist 
frequency (/at w 30000/iHz). Therefore, we can 
securely focus our analysis on the narrow region 
of the main power. The spectral region above 
20000/iHz is not plotted as it contains no addi- 
tional information. 
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Fig. 1.— Light curve of PG 0014+067 obtained 
during Xcov24. Each panel spans 24 hours, on 
the UT date indicated. Different density of data 
points comes from several factors: overlapping of 
the data, different sampling time and different 
signal-to-noise ratio. 
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WET data set. The expanded frequency range on 
the left displays the impact of the data sampling 
on the noise well beyond the 1 d~^, while the closer 
view on the right shows the details of the alias 
pattern. Both cases (top and bottom) sample a 
noise-free sinusiod; the weights for the data were 
applied to the sinusoid to compute the window 
patterns shown in the bottom panels. The vertical 
amplitude scale is amplitude (in arbitrary units). 
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Fig. 2. — Weighted Fourier transform (ampli- 
tude spectrum) of the complete WET data set on 
PG 0014-^067. 
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Fig. 4.— Weighted FT of the PG 0014-^067 WET 
data expanded to the region of the most power. 
The mean noise level (cTnoise), 3i7noise, and 3.7(7noiso 
are shown as nearly horizontal lines (from bottom 
to top respectively). 



Figure 4 expands the FT of the entire data 
set around the region of the main power, 5000 — 
10000/iHz. The sohd white hne in Figure 4 corre- 
sponds to the noise level, Cnoiso, the correspond- 
ing red (dark) line represents 3 times the noise 
level, Scnoise, while the blue (top) line represents 
3.7 times the noise level, 3.7crnoisc- The noise 
in the power spectrum is calculated by averaging 
the amplitudes around each frequency in a wide 
frequency range. For our criterion determining 
which peaks are "real" signals in the power spec- 
trum, we adopt a conservative significance thresh- 
old at 3.7(Tnoisc. Kuschnig et al. (1997) estimate 
that this results in a 99 per cent confidence limit. 
The 3crnoise threshold is adopted by others as suf- 
ficient (Brassard et al. 2001; Metcalfe et al. 2005; 
Charpinet et al. 2005) which would correspond to 
80 per cent confidence limit (Kuschnig et al. 1997). 
Hence, for comparison we show this lower thresh- 
old line (given in red) in the figures. 

Although the amplitude spectrum (Figure 2) 
does not appear very complex, attempts to de- 
termine the underlying variations by prewhiten- 
ing reveal that this star indeed has a complicated 
pulsation spectrum. After identifying the highest 
amplitude peak in the (computed) weighted power 
spectrum of the original data set, we went back to 
the unweighted FT to confirm the periodicity. We 
then removed this peak from the data by subtract- 
ing a sine wave with the frequency, amplitude and 
phase determined by non-linear least-squares fit 
(NLLS) from the light curve. In our NLLS fits, we 
use Tniax, the time of first maximum after To— 21 
August 2004 OhUT, for the phase variable. We cal- 
culate the power spectrum of the residuals (light 
curve of the entire data set with this sinusoid re- 
moved). Then we looked at this "prewhitened" 
weighted FT to find the next highest peak and 
repeated the procedure. We determined the val- 
ues of frequencies step by step, carefully checking 
whether new frequencies affected the other values 
and calculating a simultaneous n- frequency fit to 
the data until no new peaks could be identified 
with significance (with amplitude above 3.7CTnoiso)- 

Figure 5 illustrates this procedure. It shows 
the spectral window centered at 7000/iIIz in the 
top panel, and the FT of the entire data set in 
the region around the highest amplitude peaks 
in the next panel. The FT, prewhitened by the 
two highest peaks, /i — 7088.67/iHz with an am- 



plitude of 3.0 mma and /2 = 6826.06/iIIz with 
an amplitude of 2.4 mma, is presented at the 
third panel from the top of Figure 5. The high- 
est amplitude peak /i appears to have fine struc- 
ture, since significant power remains above the 
3.7(Jnoise levcl (bluc/top solid line) after remov- 
ing the highest peak at 7088.67/iIIz. Two addi- 
tional closely spaced peaks are present at frequen- 
cies /s = 7091.0AiHz and fi = 7091. 7mHz with 
amplitudes of 1.2 mma and 1.1 mma , respec- 
tively. Even though this frequency difference is 
formally unresolvable given our frequency resolu- 
tion of 0.9/;iHz, both of these modes are easily re- 
solved in our NLLS fit to the data. Furthermore, 
all the power above the 3crnoise level (red/lower 
solid line) is gone after prewhitening the data with 
those three peaks, bottom panel of Figure 5. The 
peak at 6826.06/iIIz is cleanly removed by a single 
frequency prewhitening step. 

In the extended observations of sdBV stars it 
has been detected that the individual pulsation 
frequencies have amplitudes that vary in time. 
Such frequencies could introduce spurious peaks 
in the FT as both the Fourier analysis and the 
prewhitening assume constant amplitude. There- 
fore we checked for the amplitude and the phase 
variability of the two highest amplitude periodic- 
ities (/i and /2) as we were able to resolve them 
in each individual run. The peak at 7088.67/iHz 
appears to vary both in amplitude and in phase. 
The amplitude and phase of the /2 mode does not 
appear to be changing throughout the run. The 
errors on amplitudes and phases for some runs are 
too high to allow us to make any claims. All that 
can be said is that there is a possibility of a pe- 
riodic change in amplitude and phase for the /i 
frequency with a period of roughly five days, and 
that the frequency at /2 appears to be stable both 
in amplitude and in phase over the run. Figure 6 
shows the behavior of the phase (Tmax) and am- 
plitude of /i and /2 . This figure shows the ampli- 
tude and phase of individual long runs; we grouped 
nearby short runs to reduce the phase and ampli- 
tude uncertainties when necessary. 

The variability we see in /i can be the intrin- 
sic amplitude variability of the mode or it can be 
caused by beating of two or more closely spaced 
frequencies. While the correct answer is impossi- 
ble to give within the frame of the current time- 
series techniques, all the clues of our analysis point 



toward attributing the amplitude variability in /i 
to the beating of closely spaced frequencies. If 
this was a stochastic amplitude variation, the sig- 
nature in the Fourier transform would mimic two 
or more closely spaced modes of similar amplitude, 
and complete removal with a single sinusoid in the 
prewhitening process would not be appropriate. 
Therefore, the beating of closely spaced frequen- 
cies is the most probable cause of the amplitude 
variability we see in /i as it is consistent with the 
fine structure of two additional closely spaced fre- 
quencies we found in the data. 

The next highest peak in the FT is at 7289.0/zHz. 
Figure 7 displays this region of the frequency do- 
main, showing the FT of the data before (top) and 
after (bottom) removing this periodicity, with the 
parameters determined by the NLLS fit to the 
data. Although there is some power left above 
the 3(Tnoisc level (red/lower solid line) NLLS can- 
not satisfactorily converge on the parameters of 
these low amplitude peaks. Hence, prewhitening 
stops here for this region, but we do list the high- 
est residual peak in this range as a possible real 
frequency in Table 3. 

Figure 8 shows the region of the next highest 
power; the top panel shows the residual FT after 
the five highest peaks have been removed (/i to 
/s). The highest power (amplitude) peak on the 
upper plot appears as a single peak of about 0.9 
mma. Detailed NLLS analysis finds two closely 
spaced modes at /g = 6632.8/l(Hz with an ampli- 
tude of 0.65 mma and /g = 6631.9/iHz with an 
amplitude of 0.5 mma. After those two sinusoids 
have been removed from the lightcurve, the resid- 
ual spectrum stays clean, as it is shown on the 
middle panel of Figure 8. The frequency difference 
between those two modes is of the order of our 
resolution, 0.9/iHz. The next highest peak is at 
/g = 6452.9/iHz with an amplitude just above our 
threshold (0.4 mma). Removal of this frequency 
reveals no more power above the noise (the bot- 
tom panel on Figure 8). The only power left in 
this region is at 6674/zHz, but at an amplitude 
that is below our adopted threshold, albeit above 
the 3crnoisc level. We cannot confirm this frequency 
is real, we do however indicate this frequency in 
our "almost" list (Table 3). 

At the low frequency end, there is power with a 
rather low amplitude, at about 5923^Hz. Just by 
comparing with the spectral window (Figure 3) it 
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Fig. 5. — Top panel: the spectral window centered 
at 7000/iHz given on the same frequency scale as 
the data. Second panel: the FT of the entire WET 
data set on PG 0014+067 expanded around the re- 
gion of the highest power. Third panel: the resid- 
ual FT prcwhitcned by /i and /2 at frequencies 
indicated. Bottom panel: prewhitened FT by the 
first four frequencies from the Table 2. Note the 
different amplitude scale. 
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Day (October 2004, UT) 

Fig. 6. — Top panel: phase (Tmax, in seconds) as a 
function of time for the two principal oscillations 
in PG 0014+067. Open circles are for /2 and solid 
circles are for /i. Bottom panels: amplitude (in 
mma) for /i and /2 over the course of the WET 
run. 
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Fig. 7. — Top panel: the FT of the entire data 
set, pre-whitened by /i to f^, expanded around /s. 
Bottom panel: the residual FT after prewhitcning 
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the range from 6380 to 6880 /iHz. Middle panel: 
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is clear that we are most probably dealing with 
two or more closely spaced frequencies. NLLS 
finds a mode with an amplitude of 0.5 mma at 
/t = 5923.4/zHz. Figure 9 shows this region, giv- 
ing the temporal spectrum before (top) and after 
this mode is removed by prewhitening (bottom 
panel). There is still excess power left after its 
removal, with a frequency of 5921.3//Hz. NLLS is 
unable to converge on a simultaneous fit to these 
two modes along with the nine found previously. 
This second frequency is included into the list of 
possible true periodicities (Table 3). 

The only power left above the threshold is 
just below 6200/iHz. The attempt to find the 
NLLS fit to the data including this mode was 
successful despite its low amplitude, and we have 
therefore removed this peak with a frequency of 
/io=6193.5/iHz and an amplitude of 0.4 mma from 
the total lightcurve. Prewhitening successfully re- 
moves all the excess power with only the noise left 
in that region. Figure 10 shows the FT of the 
residuals. The top panel shows the starting spec- 
trum, in this case the residual FT prewhitened by 
the nine frequencies found, and the bottom is the 
FT after we have removed the ten frequencies. 

To summarize, performing Fourier analysis and 
NLLS fitting on the entire WET data set obtained 
during Xcov24, we have detected a total often fre- 
quencies in the power spectrum of PG 0014+067 
with a 99% significance level and with a resolu- 
tion of 0.9/xHz. The final results of a simultaneous 
multiple sine-wave fit to the total light curve are 
summarized in Table 2, listing the frequencies, am- 
plitudes and phases of the main modes detected in 
the WET data set. The column labeled as Mode 
gives the name of each frequency as /„ where n 
is ordered by relative amplitude. The column la- 
beled as Tmax is the time of the first maximum af- 
ter August 21 2004 Oh UT, given in seconds. The 
uncertainties for the values given in the table are 
formal least-square estimates of the la errors. 

Figure 11 shows the FT of the combined WET 
data set on PG 0014-1-067 before (top) and after 
the removal (bottom) of ten periodicities listed in 
Table 2. We note that in the region where most 
peaks appear, the 3.7crnoisc level corresponds to 
about 0.44 mma. There are thirteen low ampli- 
tude peaks below our adopted threshold (3.7(Tnoiso) 
albeit above the somewhat less conservative accep- 
tance criterion at 3crnoisc- One of these, at 13,540 
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Fig. 9. — The residual FT of the entire data set 
at the low frequency end of the main power (top) 
after /i through /g, /s, and /g are removed. Re- 
moval of fr results in the FT shown in the bottom 
frame. 
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Table 2 
Periodicities found in WET data on PG 0014+067. 



Mode 


Frequency 


Period 


Amplitude 


T 1 




[/iHz] 


[s] 


[mma] 


[s] 


/7 


5923.4 ±0.1 


168.821 ±0.003 


0.54 ±0.13 


20 ±16 


/lO 


6193.5 ±0.2 


161.458 ±0.005 


0.44 ±0.13 


125 ±20 


/9 


6452.9 ±0.2 


154.968 ±0.005 


0.45 ±0.13 


150 ±18 


/8 


6631.9 ±0.2 


150.786 ±0.005 


0.49 ±0.15 


90 ±20 


/6 


6632.8 ±0.1 


150.766 ±0.003 


0.65 ±0.15 


143 ±15 


/2 


6826.06 ±0.03 


146.4974 ±0.0006 


2.38 ±0.13 


28 ±3 


/l 


7088.67 ±0.03 


141.0702 ±0.0006 


2.98 ±0.13 


8±3 


/3 


7091.0 ±0.1 


141.023 ±0.002 


1.22 ±0.33 


133±11 


A 


7091.7 ±0.1 


141.011 ±0.002 


1.10 ±0.32 


55±11 


h 


7289.0 ±0.1 


137.193 ±0.002 


0.65 ±0.13 


107±11 





^time of first maximum after To^ 21 August 2004, Oh UT 



fiRz, lies outside of the range displayed in Fig- 
ure 11. These peaks are at best marginal detec- 
tions. In Table 3 we list the frequencies, corre- 
sponding periods and amplitudes of those thir- 
teen marginally detected modes should they be 
detected in other data sets. 

3.2. Comparing with previous data on 
PG 0014+067 

In partial support of the run, six nights of 
observation of PG 0014+067 took place two 
months prior to Xcov24 (August 2004) using the 
high-speed multi-channel photometer ultracam 
(Dhillon & Marsh 2001) on the 4.2m WiUiam Her- 
schel Telescope (WHT) at La Palma (Jeffery et al. 
2005). We compare our results with published re- 
sults on PG 0014+067 from Brassard et al. (2001) 
and Jeffery et al. (2005). Table 4 hsts the pe- 
riodicities identified by Jeffery et al. (2005) and 
Brassard et al. (2001) together with the WET 
results. 

Note that all frequencies identified by the WET 
(Table 2) were identified by either Jeffery et al. 
(2005) or by Brassard et al. (2001) or by both 
to within the cycle per day aliases. The high- 
est amplitude frequency in the WET data is at 
7088.67/iHz and it is seen as a multiplet of three 
closely spaced frequencies (along with 7091.0/iHz 
and 7091.7^Hz). The ultracam data set identi- 
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Fig. 11.— Total FT of the original WET data 
set on PG 0014+067 (top) followed by the resid- 
ual FT of the same data set, prewhitened by the 
ten frequencies listed in Table 2. Note that the 
red (lower) solid line corresponds to 3crnoisc, while 
the blue (upper) solid line corresponds to 3.7crnoise. 
The white solid line is Cnoiso. 
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Table 3 
Marginal detections- peaks with an amplitude between 3.7cr„ 

DATA ON PG 0014+067. 



Frequency 


Period 


Amplitude 


[mHz] 


[s] 


[mma] 


5921.3 


168.882 


0.40 


6674.0 


149.835 


0.43 


7278.9 


137.383 


0.40 


7895.0 


126.662 


0.39 


7938.3 


125.971 


0.45 


8870.8 


112.729 


0.42 


9576.6 


104.421 


0.41 


9834.2 


101.686 


0.38 


10807.8 


92.526 


0.36 


11116.5 


89.956 


0.36 


11364.8 


87.991 


0.36 


12083.7 


82.756 


0.37 


13539.8 


73.856 


0.38 



fied 4 closely spaced frequencies in this range (see 
Table 4) out of which the highest amplitude one 
in their data set is at 7091. 7^Hz. This multiplet 
is not identified in the CFHT data probably due 
to the resolution of their run, only one frequency 
of this multiplet is found at 7088.7/iHz. The am- 
plitude and phase of the highest amplitude mode 
(7088.67/xHz ) in the WET data set as well as the 
highest amplitude mode (7091.7/iHz) in the UL- 
TRACAM data set were noticed to vary on the same 
time scale of five days. 

The amplitudes of the individual modes in 
sdBV stars are known to change on the time- 
scale of years, months and even days (Kilkenny 
et al. 1999), although the question whether this 
results from true physical change in the star or 
beating between closely spaced frequencies still 
remains unanswered. In this particular case both 
the WET and ultracam data sets find the ampli- 
tude and phase variability on the same time scale 
and identify closely spaced frequencies, a multi- 
plet around 7090/iHz. This gives an additional 
credibility to the hypothesis that this relative am- 
plitude change is caused by beating, though longer 
runs would be needed to confirm this speculation. 



4. Discussion 

The frequencies we saw in PG 0014+067 during 
the October 2004 WET campaign showed no obvi- 
ous patterns that could readily be attributed to ro- 
tational splitting or other common effects in non- 
radially pulsating stars. However, some unusual 
systematics among the frequencies became read- 
ily apparent. In this section we discuss those sys- 
tematics, and explore the statistical significance of 
such patterns by including additional secure fre- 
quencies seen by Jeffery et al. (2005). 

4.1. Strange systematics in the observed 
frequencies 

The frequency list of PG 0014+067 displays 
some suggestive systematics. As an exercise in 
numerology, but one that will eventually provide 
some statistically significant results, we examined 
these systematics further. 

Considering only the dominant modes in each 
closely-spaced multiplet, there are 7 main frequen- 
cies in the WET data on PG 0014+067, repre- 
sented by /i, /2, /s, /e, /7, /g, and /iq. Within 
this group, we find several pairs of frequencies sep- 
arated by multiples of 90/xHz. The first five lines 
in Table 5 list these differences. While not all en- 
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Table 4 

WET FREQUENCIES IN PG 0014+067 (in ^Hz) compared with the ultracam data, Jeffery et 

AL. (2005) and CFHT data, Brassard et al. (2001). 



Mode 


WET 


ULTRACAM 


CFHT 






5780.9" 


5896.2 


/7 


5923.4 


5921.9" 
5924.8 


5923.2 


/lO 


6193.5 




6227.7 = /lo + 3d-^ 


/9 


6452.9 


6454.4 




/8 


6631.9 




6630.7 


/6 


6632.8 


6632.6 
6646.5 
6659.9" 
6726.8" 


6621.1 = /s - Id-i 


/2 


6826.06 


6826.1 


6837.5 = /2 + ld-^ 






7076.6" = /i-ld-^ 


7079.1= /3 - Id "^ 


h 


7088.67 


7089.1" 


7088.7 


h 


7091.0 






h 


7091.7 


7091.67* 
7093.4 
7094.8 

7187.5 


7150.2 


h 


7289.0 


8588.9" 


7286.2 
7670.3 
7952.1 
8552.1 

9797.6 






9971.5 


9970.3 






11547.9 




(2 X /lo) 






12386.8 


(2 X /g)" 




12910.9" 





* highest amphtudc seen in ULTRACAM data 

^secn in white light only 

"^niultiplct 

^/q as seen in ULTRACAM data 
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tries in the table are independent, it illustrates 
that some (large) frequency differences are very 
similar to one another, and all share a common 
factor of 90.47 ±0.53 mHz. 

The frequency f^ does not pair with any other 
mode in a difference that is a factor of 90.47 
/zHz. However, the average value of f^ and /i, 
7188.8/iHz, is four times that interval from /2. 
This somewhat suspicious additional entry in Ta- 
ble 5 is supported by the observation of a fre- 
quency that lies within 1.3/^Hz of that value by 
Jeffery et al. (2005) (the frequency of 7187.5/iHz 
in Table 4). Thus the three frequencies {f^ = 
7289.0AiHz, 7187.5AiHz, and /i = 7088.7/iHz) form 
a nearly equally-spaced triplet with a spacing of 
about 101.5/iHz. In light of this, we explored the 
frequencies in Table 2 allowing for an « 101/iHz 
deviation; these are shown in the last four lines of 
Table 5. With this additional complication, addi- 
tional evidence for a possible chain of frequencies 
separated by about 90/iHz appears. If we adopt 
a value of 90.47/iHz for the first spacing, then the 
last four lines of Table 5 suggest that the second 
spacing is 101.1±1.4 /iHz. 

We recognize that the spacings identified by 
this subjective technique alone may not be unique 
or even meaningful. We note that there is no a 
priori expectation of a chain of equally spaced (in 
frequency) modes for these stars, nor is there an 
expectation, based on the physics of the star, for 
a splitting as large as lOl/iHz. Neither 90/xIIz nor 
100/iHz spacings are seen in theoretical models 
of these stars. The asymptotic frequency spac- 
ing for p— modes in these stars (that is, modes 
with n > I, consecutive n, and alternating I) is 
roughly 750/iIIz. This value was determined using 
the model frequencies from Brassard et al. (2001) 
along with our own sdB models with appropriate 
values for logg and Tcff- 

Even so, could the lOl/iHz splitting be a rota- 
tional splitting? While this large splitting could 
be caused by rotation, the implied rotation rate 
(if solid-body rotation) although not ruled out, 
would be faster than any known single sdBV de- 
tected so far, based on spectroscopic study of fine 
profiles. Spacings this large are seen only in close 
binary pulsating sdBV stars such as PG 1336-018 
(Kilkenny 2003). On the other hand Kawaler & 
Hostler (2005) suggest that rapid internal rota- 
tion could produce large splittings in a star with 



a slow surface rotation rate, but their predictions 
suggest that the splittings of modes differing in n 
and / should not show the same value. 

4.2. An empirical, phenomenological rela- 
tion 

With two apparent splittings present, we de- 
cided to explore an entirely phenomenological pa- 
rameterization that could then be used to make 
an empirical fit to the observed frequencies, and 
provide a framework for determining the statisti- 
cal significance of these patterns. We chose a form 
reminiscent of asymptotic p— mode pulsation with 
a constant rotation frequency 



f{i,j) = fo + i^S + jxA 



(1) 



where 6 represents a small spacing (and i is an in- 
teger ranging from upwards) and A represents 
a large spacing (with j initially limited to being 
either -1, 0, or 1). In the equation above, fo rep- 
resents a zero-point for the fit with i = j = 0. 

In the general case of fitting an observed set of 
frequencies, we performed a two-dimensional x^ 
minimization over a grid of values of S from 40 to 
120 fiRz and A from 80 to 210 fiRz. For each {6, 
A) pair, we found the combination of {i, j) that 
minimized the root-mean-square (RMS) difference 
between the observed and model frequencies, and 
computed the x^ of the fit. Clearly, there is an 
aliasing problem when the combinations of i, j, 
S, and A produce commensurate spacings, so the 
X^ surface shows multiple minima. To help break 
that degeneracy, we impose an additional crite- 
rion that a ((5, A) solution result in at least two 
modes that have the same value of i but different 
j. The quantized nature of the problem (the val- 
ues of i and j are constrained to be integers) and 
the additional constraint on the values of j com- 
plicate determination of statistical significance of 
any fits. Therefore to find the form of the distribu- 
tion of x^ for a random selection of frequencies, we 
performed a series of Monte Carlo trials using fre- 
quency lists drawn from uniformly random distri- 
butions within the range shown by PG 0014+067. 

Using only the WET frequencies, we find two 
solutions of high significance; however, only one 
satisfies the requirement of having two modes with 
the same value of i but different values of j . That 
solution is the one found by the successive differ- 
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Table 5 
Frequency differences in the PG 0014+067 WET data 



Frequency pair 


difference 




[mHz] 


/6 - /9 


179.9 


/lO - /t 


270.1 


/i " /e 


455.9 


/2 - /lO 


632.6 


/i - /g 


635.8 


(/5,/l)-/2 


362.8 



multiple 

[mHz] 



/5 - /l 

/l - /2 
/s — /2 

/q — ./lO 
/2 - /9 



200.3 

262.6 
462.9 
259.4 
373.2 



== 2 X 89.95 
== 3 X 90.03 
= 5 X 91.17 
= 7 X 90.37 
= 7 X 90.82 
= 4 X 90.69 

== 2 X 100.15 



: 4 X 90.47^ 
4 X 90.47^ 
4 X 90.47^ 
3 X 90.47^ 



- 99.29 
+ 101.01 

- 102.49 
+ 101.78 



mean spacing — 90.47+ 0.53 



mean spacing = 101.1+1.4 



^fixed at mean value from lines 1-5 



encing procedure described in the previous sub- 
section. The fit values were S — 90.37 /iHz, 
A = 101.46 iMz, and /o=5922.46 iiRz. 

With only seven frequencies and more than 
three free parameters (the values of 6, A and fo are 
free, but the choices of i and j are constrained), 
the fit above is not well constrained. Via the 
Monte Carlo technique with 7 frequencies, we find 
the fit to be significant at the 95% confidence level; 
that is, in a trial of 285 frequency sets, only 13 had 
a lower value of x^ ^^d at least two modes with 
the same value of i and different values of j . While 
this is encouraging in terms of suspecting that the 
parameterization above may be meaningful, it is 
not, by itself, very convincing. 

Fortunately, we have additional frequencies in 
this star that have been reliably measured by oth- 
ers that lie within the frequency range from the 
WET data alone. This list is given in Table 6. 
By "reliable" we mean peaks seen in more than 
one investigation, and those that are clearly true 
peaks and not aliases. Given the exceptional qual- 
ity and clean window of the ultracam data, we 
give that data set higher relative weight, using the 
WET to resolve residual ambiguities originating 



from diurnal aliases. Using these additional ob- 
servations, five additional frequencies in the range 
that the WET showed pulsations can be examined 
for compliance with the relation derived using only 
the seven WET frequencies in Table 4. 

Of the five additional frequencies in Table 6, 
four fit the pattern determined from the WET fre- 
quencies alone using the same parameters. Again, 
we use Monte Carlo simulations of this process 
to judge whether or not this is a significant ef- 
fect. In this case we measure the RMS deviation of 
five additional randomly chosen frequencies from 
the previously-established relation from the WET 
data. The RMS deviation is calculated after drop- 
ping the worst-fitting of the five. Comparing the 
distribution of the RMS deviations in the Monte 
Carlo trials with that from the data in Table 6 re- 
veals that the fit of the relation (equation 1) with 
the data is significant at the 99.4% confidence level 
(based on 9829 trials) under these conditions. 

In practice, we impose an additional require- 
ment that the fit must include modes that have at 
least two pairs with the same i but different j (i.e. 
that at least two values of i show more than one 
j component) . Enforcing this additional condition 
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on the fit criteria restrict the ahasing problem and 
result in a smaller number of possible solutions. 
This process also results in a much higher statis- 
tical significance for the fit to PG 0014+067; the 
confidence level rises to the 99.95% level. In the 
fit in Table 6, we see that there are three values of 
i that have more than one j component, with one 
complete 'triplet' - a concordance that none of the 
9829 Monte Carlo trials reached. 

Finally, we can employ the fit procedure us- 
ing all of the modes in Table 6 to determine more 
accurate values of the parameters of Equation 1. 
For PG 0014-1-067 this procedure yields values of 
(5=90.37, A=101.22, and /o=5923.24 (aU in /xHz). 
Using these values in Equation 1, we find the 
model frequencies listed in Table 6. Note that the 
fit is extremely good - all but one of the 12 modes 
identified in the table are fit to within 0.05%, with 
an RMS difference of 0.013 5. To illustrate the 
closeness of the fit. Figure 12 shows the frequen- 
cies of PG 0014+067 in an "echelle diagram" sim- 
ilar to those used in helioseismology. The vertical 
axis is frequency, and the horizontal axis (also in 
frequency units) shows the departure of each mode 
from an integral multiple of the small spacing 5. 
The repeating pattern of j = modes stacks above 
an ordinate of 0, while modes with nonzero values 
of j lie on either side at frequencies A away. 

Is this asymptotic pulsation? After all, high- 
order p-modes show a more-or-less constant fre- 
quency spacing; such sequences are seen in he- 
lioseismic data and in the rapidly oscillating Ap 
stars (i.e. Kurtz et al. (2005)). The sequence 
of modes split by integral multiples of 5 cannot 
be asymptotic p— mode behavior. As mentioned 
above, models of PG 0014+067, and sdBV pul- 
sators in general, indicate that the radial funda- 
mental frequency in the models is usually close 
(in frequency) to the observed mode frequencies. 
Asymptotic relations such as Equation (1) are usu- 
ally valid (at the few-percent level) only for values 
of n ;3> ^ or more generally for large values of n. 
Even so, the computed frequency separation for 
p— modes in sdBV models yields values of several 
hundred /xHz - a factor of 5 or more larger than 
what PG 0014+067 shows. 

In summary, the frequencies seen in PG 0014+067 
obey a simple empirical trend - one that cannot 
be understood in terms of standard evolutionary 
and pulsation models for such stars. Simply put. 



there is no applicable physics in this purely em- 
pirical fit. In a subsequent paper, we will examine 
the frequencies present in other well-studied sdB 
stars, and we show that they show a similar pat- 
tern as well. The theoretical interpretation of this 
pattern is not at all clear, but we anticipate that 
the systematic frequency differences discovered in 
PG 0014+067 and other sdB stars could be a very 
important clue to unlocking some secrets of their 
interiors. 
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Table 6 
Pulsation frequencies in PG 0014+067, along with an empirical model fit. 



# 


Frequency 


Amplitude 


note 


i 


J 


Model 


Difference 




[mHz] 


[mma] 








[mHz] 


[a*Hz] 




5780.9 


0.35 


ULTRACAM (w) 










/7 


5923.4 


0.54 


fine structure 








5923.2 


0.2 


/lO 


6193.5 


0.44 




3 





6194.1 


-0.6 


/9 


6452.9 


0.45 




7 


-1 


6454.1 


-1.2 


/6 (/s) 


6632.8 


0.65 


fine structure 


9 


-1 


6634.6 


-1.8 




6646.5 


0.60 


ULTRACAM 


8 





6645.6 


0.9 




6659.9 


0.34 


ULTRACAM (w) 


7 


1 


6656.5 


-3.4 




6726.8 


0.37 


ULTRACAM (w) 


10 


-1 


6724.9 


1.9 


/2 


6826.1 


2.38 




10 





6826.1 


0.0 


/l (/3,/4) 


7088.7 


2.98 


fine structure 


14 


-1 


7086.1 


2.6 


- 


7187.5 


0.66 


ULTRACAM 


14 





7187.3 


0.2 


/5 


7289.0 


0.65 




14 
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Fig. 12. — "Echelle" diagram of the frequencies in 
PG 0014-1-067; the frequencies have been folded on 
the small spacing 5 = 90.37 and stacked, showing 
the uniformity of the frequency spacings within 
the main band and the two side bands separated 
from the main band by the large spacing A — 
101.22. Expanded points are from the non-WET 
modes in Table 6. 
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